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Abstract. We consider the application of the magnetic flux leakage (MEL) method to the detection 
of defects in ferromagnetic (steel) tubulars. The problem setup corresponds to the cases where the 
distance from the casing and the point where the magnetic field is measured is small compared 
to the curvature radius of the undamaged casing and the scale of inhomogeneity of the magnetic 
field in the defect-free case. Mathematically this corresponds to the planar ferromagnetic layer in a 
uniform magnetic field oriented along this layer. Defects in the layer surface result in a strong de¬ 
formation of the magnetic field, which provides opportunities for the reconstruction of the surface 
profile from measurements of the magnetic field. We deal with large-scale defects whose depth is 
small compared to their longitudinal sizes—these being typical of corrosive damage. Within the 
framework of large-scale approximation, analytical relations between the casing thickness profile 
and the measured magnetic field can be derived. 
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1. Introduction 

The magnetic flux leakage (MFL) method is a powerful tool for non-distructive inspection of the 
integrity of ferromagnetic casings 03 El 12 01 [5. 61 0 or, more generally, for determining the shape 
of ferromagnetic objects. The basic idea of the method is to reconstruct the shape features of the 
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Figure 1: Mirror-symmetric ferromagnetic layer of non-uniform thickness in the magnetic field 
H and the coordinate frame 

inspected ferromagnetic object from the deformation of the magnetic field. The conventional way 
of applying the MFL method is (i) recognition of typical patterns in the magnetic field measure¬ 
ment data, (ii) identification of typical defects of the casing corresponding to these patterns, and 
(iii) evaluation of geometrical parameters of the defects identified from the patterns. Such an ap¬ 
proach faces obvious problems when one deals with complex-shape defects such as metal losses 
due to corrosion. Although the problem of the reconstruction of an arbitrary object shape from 
measurements of the magnetic field (all three components of the H -field should be measured) on 
some surface in space near the object is mathematically well-posed and resolvable (numerically), 
in practice one encounters issues which make this approach generally inapplicable (which is the 
reason for using the conventional practice outlined above). However, specifically in the case where 
complex-shape defects are most important—in the case of corrosion damage—there is an oppor¬ 
tunity for reconstruction of the arbitrary casing thickness profile based on an analytical approach. 
This analytical approach is possible owing to the model reduction—the large-scale approximation, 
which assumes the defect depth and the casing thickness to be small compared to the defect width, 
as it is typical for corrosion metal losses. The corrosion loss of metal can be non-large-scale near 
a weld, where corrosion rapidly advances along the contact interface. Otherwise, the large-scale 
approximation is reasonable. 

In this paper we consider the application of the MFL method to the detection of defects in 
ferromagnetic (e.g. steel) tubulars including wellbore casings. Our treatment is focused on the ac¬ 
curate reconstruction of casing thickness profiles, in contrast to the conventional approach which 
is the recognition of magnetic field patterns corresponding to catalogued typical defects by means 
of neural networks or similar data analysis tools. The problem setup we use corresponds to mea¬ 
surements made with modem devices designed for the MFL inspection of wellbore casings (e.g., 
see 0). In this setup the distance from the casing and the point where the magnetic field is 
measured is small compared to the curvature radius of the undamaged casing and the scale of in¬ 
homogeneity of the magnetic field in the defect-free case. Mathematically, this corresponds to the 
planar ferromagnetic layer in a uniform magnetic field oriented along this layer. Defects of the 
layer surface result in a strong deformation of the magnetic field, which provides opportunities for 
the reconstruction of the surface profile from the measurements of the magnetic field. 
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2. Analytical theory: Large-scale approximation 

We restrict the analytical treatment to the case of large-scale defects, i.e. defect length and width 
is large relative to the depth or the casing (metal layer) thickness. For this case we show that the 
magnetic field is sensitive to the casing thickness profile, but not to the inner and outer surface 
profiles independently. Hence, as a starting point, one can address the problem of a symmetric 
layer, one surface of which is a mirror image of the other. 

2.1. Mirror-symmetric layer 

Mathematical description of the problem 

Let us consider the ferromagnetic layer confined between two surfaces z = £(x, y) and z = 
— ((x, y ), where the (x, y) -plane is the middle plane of the layer and the z-axis is orthogonal to 
it. The uniform external magnetic field H is applied along the x-axis. The system is sketched in 
Figure Q] 

We adopt the following assumptions for the problem: 

1. The layer geometry and fields possess the symmetry property (z —* —z); 

2. The linear magnetisation law for both the ferromagnetic material and the material around it 
is given by: Bj = HjHj, j = 1,2; 

3. The magnetic permeability of the surrounding material is small compared to that of the 

ferromagnet, — 1; 

F2 

4. Surface defects are large-scale, which means the typical longitudinal size of defects L > ( 
and, therefore, |V 2 CI ^ 1- 

(In the following subsections we will extend our consideration beyond restrictions (1) and (2).) 
Henceforth, the subscripts of fields and parameters indicate the corresponding domain (1: ferro¬ 
magnet, 2: upper outer area, 3: lower outer area); for the gradient and Laplace operators, V 2 and 
A 2 , the index “2” indicates the two-dimensional versions of them calculated with respect to x and 
y coordinates only. The ranges of parameter values of practical interest are presented in Table Q] 
and are consistent with the assumptions made. 

According to MaxwelFs equations, we have the following equation system 

j V x Hj = 0, 

l v-4 = o, 

with boundary conditions 

H It Ff 2r . B\ n f? 2 n 


3 



A.V. Pimenova, D.S. Goldobin, et al. 


Magnetic Flux Leakage: Analytical Theory 


Table 1: Reference values of real system parameters 


ratio of magnetic permeabilities 
longitudinal scale of defects 
location of magnetic sensors 


/11//12 : 100 - 1000 

L : (10-20)C 
2 : (2 - 8)C 


for the normal to the surface (subscript “n”) and tangential (“r”) components of magnetic field, 
respectively. When the curl of a vector field is zero in a simply-connected domain, one can in¬ 
troduce the scalar potential for this field within this domain. We introduce the scalar potential $, 
H = — V<T, which obeys the equation 

A$ J = 0, (2.1) 

while the boundary conditions read 


lii 


$1 = 

dn 


$ 2 




d § 2 

^ 2 dn 


( 2 . 2 ) 

(2.3) 


Since the magnetic permeability of the ferromagnet is considerably larger than that of the 
surrounding material, to the leading order of approximation the flux of the magnetic field does 
not go outside the boundaries of the ferromagnet, i.e. 


C 


dx 


dz I = 0 . 


(2.4) 


In this case the normal derivatives of <3> on the boundary are equal to zero, and Eq. (12.31) takes the 
form 

£ = 0. (2.5) 

dn 

For infinitely large scale inhomogeneities the magnetic //-field within the layer is the same as 
for the defect-free planar layer, Hi = H. Hence, we can look for the correction to the uniform 
field H. One can write down the Taylor series for $1 with respect to z 

2 4 

$ 1 0, y, z ) = -Hx + $i 0) (x, y) + $i 2) (x, y)^ + {x,y)^ + ... 


(only even powers of z are present due to the symmetry z —> —z). Substituting this series into 
Eq. (12.11) and renaming T, (x, y) as F(x, y), one finds 


$i(x, y, z) = -Hx + F(x, y) - A 2 F(x, y)^ + A\F(x, y) '— - ... . (2.6) 
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Notice, here the z 2n -term is of the order of magnitude of / L) 2 " and thus only first several terms 

can be important for the large-scale case. Since F vanishes for infinitely large scale of defects, it 
should be small compared to the leading term of $1 for finite large scale L by continuity, i.e. 


|v 2 f| < n. 

On the surface 2 = ((x, y), Eq. (12.21) yields 


$ 2 (* = 0 


$1 = — Hx + F 


-A 2 FC 2 + o 



(2.7) 


2.2. Two-dimensional case 


Let us consider the two-dimensional problem of a ferromagnetic layer uniform in the ^/-direction. 
One can see that for the two-dimensional case the integral in Eq. (12.41) can be found as 


1: 


(M> 1 
dx 


dz = const = 


( 2 . 8 ) 


where Co is the ^-coordinate of the undamaged surface. This significantly simplifies the task and 
makes it possible to solve the problem analytically. Substituting expression (12.61) into the latter 
equation, one can see that 


<I> 2 can be expanded into a series near the surface 


Cl 

' L 2 


<h 2 (C) = ®2{h) + 




dz 


(C - h) + 


1 


z=h 


2 dz 2 


0 C~h) 2 + 


(2.9) 


( 2 . 10 ) 


z=h 


To calculate (d 2 ^ 2 /dz 2 ) one can employ Eq. (12.11) . 

<9 2 <f> 2 d 2 $n 


dz 2 


z=h 


dx 2 


( 2 . 11 ) 


z=h 


Hence, Eq. (12.101) can be rewritten in the form 


$ 2 (h) = -Ux + F- 


1 d 2 F 

2 dx 2 


c 2 + 


d&o 


dz 


(h~ 0- 


1 a 2 <f> 2 


z=h 


2 dx 2 


(h-0 2 + O (F 


z=h 


Xh-CY 


L 3 


( 2 . 12 ) 

Substituting Eq. (12.91) and differentiating the last equation with respect to x, one can evaluate the 
cc-component of the magnetic //-field measured at the height h above the layer; 


H\ -7/ Co I 9 f 9 ® 2 
xlz=h ~ H C dx { dz 


z=h 


(C - ft)) + o (n^j 


(2.13) 
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Let us seek a series expansion for ( = Ci + C2 + • • • > where each term of the series is small 
compared to the previous one. To the leading order of accuracy the last equation yields 


Cl 


Co 


n 

H x \ z=h 


Substituting (12.141) into Eq. (12.131) and considering H = —V$, one can obtain 


C2 


UQ d 

H l\z=h dx 



Finally, to the terms of order ('HQ 2 / L 2 ) 


c 


Co 


n 


1 


Co d 

H x \ z=h dx 



(2.14) 


(2.15) 


(2.16) 


2.3. Asymmetric ferromagnetic layer 


In this subsection we argue that for large-scale defects the case of an asymmetric layer is equivalent 
to the case of mirror-symmetric layer. Let us consider the layer confined between surfaces z 2 = 
C2OE, y) and ;j 3 = —QQx, y ) with ( 2 Q (3- It I s convenient to introduce the middle surface z rn = 

Cm(x,y). 

> _ C2 — C3 

U 2 ’ 

and use the coordinate frame 

X = X, Z = Z-( m {x). 


Then 


d d d( m d 


dx dz dx dz 

In the new coordinates Eq. (12.11) takes the form 


and 


d_ _ d_ 
dz dz ' 


A<h 


A $i- 


9 2 Cm 2 ^Cm_^_^ , ( ^Cm \ 
dx 2 dx dx) \ dx ) 


d 2 $i 

di 2 ' 


(2.17) 


For large-scale defects the additional terms of Eq. (12.171) affects Eqs. (12.41) and (12.101) only in high- 
order terms. Hence, Eq. (12.131) is not affected by asymmetry. 


2.4. Nonlinear magnetisation law 

Ferromagnetic materials are characterized by hysteresis and non-linearity of the magnetisation. 
Hysteresis is to be avoided in inspection applications as it leads to a loss of uniqueness in the 
solution of the profile-reconstruction problem. For this reason, strong constant magnets are used 
in practice, and the system operates under conditions close to magnetic saturation (see Figure [2]). 
Close to saturation hysteresis becomes insignificant whereas the non-linearity becomes pronounced 
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Figure 2: Schematic dependence of the magnetic field B on the magnetic //-field in ferromagnetic 
material (black solid curves) and its linearization in the working parameter range (red line) 


and influences applications of the MFL method 0. Nonetheless, non-linear-magnetisation prob¬ 
lems can be treated analytically for large-scale defects as well, because the magnetic //-field within 
the ferromagnetic layer deviates slightly from H. For brevity, we consider the case of a mirror- 
symmetric layer in this section. 

As V x H = 0, we still can use substitution H = — V<f>. For B = the equation 

V • B = 0 takes the form 


V • // + 


-L°£.a.v* 

fiH OH 2 


0 . 


Substituting the series 


<f>i(x, z) = —Bx 




+ $i 2) (x)^ + $i 4) (x)^ 


+ ... 


into the latter equation and collecting z-free terms, one finds 

a 2 ®? 1 . t (2> ld^ dH* 

x 2 1 /ii// dH dx 2 

With H 2 = (H — (d / dx)<&i ]) +... ) 2 = H 2 — 2'H{d / dx)^^ +...; to the leading order of accuracy, 
one can obtain 


$j 2) = 
P = 


— (1 + / 3 ) 


<9 2 <f>i 0) 
dx 2 


+ ... , 


H d/j,i 
/ii {%) ~dH 


H=n 


(2.18) 

(2.19) 


One can show that 0 < (3 < 1. Indeed, let us write B = a(H + H 0 ) for H next to B (see 
Figure©. Then /3 = H 0 / (V, + H 0 ) with positive // 0 ; therefore, d G (0,1). On the hysteresis loop, 
which is out of the scope of our study, /? can be beyond this range. 
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Eq. (12.21) holds valid for nonlinear magnetisation and yields a modified version of Eq. (12.61) : 

1 n2 (f)(0) 

<E 2 (z = C) = $1 = -Hx + $! 0) - - (1 +p) -Q^rC + • • • • (2.20) 

Since 0 < /? < 1 for real ferromagnetic materials, the change of a coefficient ahead of the 
last term does not change its order. As was demonstrated for the case of a mirror-symmetric 
layer, any terms of this order do not affect the solution to the concerned accuracy. Therefore the 
nonlinearity of magnetisation does not affect the leading order of the equations we suggest for the 
profile reconstruction procedure for large-scale defects. 


2.5. Three-dimensional ferromagnetic layer 

For the case of the three-dimensional layer Eq. (12.41) can be solved only up to the gradient of an 
arbitrary harmonic function of x and y, which is not very helpful for our purposes; the problem 
requires a somewhat different approach compared to the two-dimensional case. 

Let us denote the magnetic potential at the height £ 0 as $ 2 0 • One can write down the Taylor 
series for $ 2 at the height 

(C - Co) + 2 -qJ- (C - Co) 2 + • • • • (2.21) 

z=Co 2=Co 

Taking into account the boundary condition <3>! = <f> 2 , one can equate Eq. (12.61) with the latter 
equation and evaluate the leading-order term of the magnetic potential within the layer: 

(C -Co ) + o(n^y (2.22) 

Z=Co 

Substituting this series into Eq. (12.41) . one finds 

V 2 • (c H + C(C - Co)V 2 /L + CV 2 C H z ) = 0 , (2.23) 

where H = — V4> 2 | z= ^ 0 is the magnetic field at the height £ 0 and H z is its ^-component. To 
solve this equation numerically it is convenient to use the exponential representation of ( = 
( 0 e~ a(0) ~ a(1> ~'". To the first two orders Eq. (12.231) takes the form 

- CV 2 a(°) • H - CV 2 aW . H + CV 2 ■ H + V 2 [C(C - Co)V 2 iT a + (V 2 (H Z ] = 0 . (2.24) 


$i(x,y,z) = $ 2 | z=?0 + 




dz 


$ 


2|- =C - ®2 


+ 




dz 


Hence, 

V 2 a(°) -H = X7 2 -H, (2.25) 


V 2( xW • H = Co[(l - 3e-^)X7 2 aW . V 2 H 2 - (1 - e~ a(0) )A 2 H z 

+ 2 e~ aW {V 2 a (0) ) 2 H z - e - a(0) A 2 a^H z \ . 


(2.26) 


Notice, <7® is not necessarily small. 

With Eqs. (12.251) and (12.261) one can calculate the layer thickness profile = 2( 0 exp(—cr (0) — 
— ...) from the magnetic field H (or some of its components) measured at non-large elevation 
above the layer. 
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Figure 3: (a): The reconstructed profile (12.161) of a ferromagnetic layer (red dashed line) com¬ 

pared to the original profile (black solid line) for the parameter values Co = 0.1, a = 0.01, 
k = 2-7T, //i / fj '2 = 100. (b): Inaccuracy (the deviation from the original profile) of the profiles 
reconstructed from the solution with direct numerical simulation (DNS) and the linear-in-defect 
solution in Fourier space. 

3. Application of the analytical technique 

3.1. Validation of applicability of the analytical technique with results of 
numerical simulation for two-dimensional case 

In order to validate the applicability of the analytical results derived, we have considered the model 
case of a ferromagnetic layer of /r 2 /fti = 100 with profile C = Co + « cos kx with k = 2n, Co = 0.1, 
a = 0.01. The magnetic field for this case was calculated both 

• with direct numerical simulation, employing a finite volume method and the mesh size dx = 
dz = 0.01, and 

• analytically in Fourier space within the framework of the linear-in-defect approximation. 
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Linear in (£ — ( 0 ) solution to the problem can be found analytically in Fourier space. Non¬ 
diverging for z —* oo solution to the problem reads 


<f)(x, z) = 


—Kx + A 


cosh k\z\ 


sin kx, \z\ < C ; 


cosh k (o 

—Kx + Ae~ k ^~^ sin kx, \z\ > ( . 


(3.1) 


Here 


A = 


(/ii - H 2 )Ka 


Hi tanh k ( 0 + /x 2 e 

One can differentiate the latter solution to find the components of the magnetic field at height h\ 


H x \ z= h — K — kA e ki ' h ^ cos kx, 
H z \ z= h = kA e~ k ( h ~^ sin kx . 


(3.2) 


In Figure [3j with synthetic data from direct numerical simulation, one can see the surface 
profile can be well reconstructed with Eq. (12.161) . The accuracy of analytical solution (13.11) and the 
role of correction (compare Eqs. (12.141) . (12.151) and (12.161) 1 can be judged from Figure[3lb). 


—> 

3.2. 3D layer: measurable //-field 

For the case of sensors measuring the H- field at certain elevation above the layer with a dense 
enough grid of measurements points (e.g., 0), one can directly employ Eqs. (12.251) and (12.261) 
with //-field derivatives approximated by finite differences. If the elevation height is of the same 
order of magnitude as the characteristic defect width, one needs first to calculate //-field on the 
“imaginary” undamaged surface and then use this field for calculation of the layer thickness profile. 
For this calculation one can use the Taylor expansion of //-field. Indeed, all the ^-derivatives of 
//-field can be calculated from x- and ^-derivatives of measured fields, because H is a gradient of 
a harmonic function (see Eq. (12.11) 1: 

d H z 
dz 

d 2 H z 
dz 2 

d n+2 H z 
dz n+2 

0H X 
dz 

d n+2 H x 
dz n+2 


dH x dH„ 


dx dy 

d 2 d 2 

+ 


dx 2 dy 2 


H z , 


d 2 d 2 \ d n H, 

+ 


dx 2 dy 2 ) dz 


n = 1, 2, 3,.. 


dH z 

dx 


d 2 d 2 \ d n H 

+ 


dx 2 dy 2 ) dz n 


n — 0,1,2, 
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dHy _ OIL 


dz dy 


d n+2 H y f d 2 + d 2 \ d n H y 


dz n+2 \dx 2 dy 2 J dz n 


n = 0,1, 2,... . 


3.3. 3D layer: measurable H z or ( dH z /dx ) 

When only one component of the magnetic field is measured, one has to make a more substantial 
use of the harmonic property of the field. Indeed, a harmonic function within 0 < x < L x and 
0 < y < L y , bounded at z — y +oo, can be represented in the basis of exponentials; 



(3-3) 


Amplitudes (j) mn can be evaluated from Fourier decomposition of the measured component of the 
H- field (or its x-derivative) 



(3.4) 


With (fimn known, one can calculate derivatives of <f>(x, y. z) at any z, which are components of the 
//-field, and employ Eqs. (12.251) and (12.261) . 

4. Conclusion 

We have developed a technique for the analytical calculation of ferromagnetic casing thickness 
profiles from measurements of the magnetic field above the layer when a homogeneous exter¬ 
nal magnetic field is applied, i.e., for the magnetic flux leakage (MFL) method for inspection of 
wellbore casing integrity. The analytical results have been derived within the framework of the 
large-scale approximation of defects, the widths of which is large compared to their depth and 
layer thickness; this approximation is generally relevant for corrosion damage (with the exception 
of corrosive damage of welds). The technique has been shown to be applicable for a nonlinear 
magnetisation law and without hysteresis within the working range of //-field strength. The lat¬ 
ter restriction potentially diminishes the applicability of the result, but MFL tools are designed to 
saturate the casing to minimize the impact of hysteresis on the analysis. The applicability of the 
analytical results has been validated with the results of direct numerical simulation. 
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